Reconstruction of a scalar-tensor theory of gravity 
in an accelerating universe 
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The present acceleration of the Universe strongly indicated by recent observational data can be 
modeled in the scope of a scalar-tensor theory of gravity. We show that it is possible to determine 
the structure of this theory (the scalar field potential and the functional form of the scalar-gravity 
coupling) along with the present density of dustlike matter from the following two observable cosmo- 
logical functions; the luminosity distance and the linear density perturbation in the dustlike matter 
component as functions of redshift. Explicit results are presented in the first order in the small 
inverse Brans-Dicke parameter lu~^ . 



PACS numbers: 98.80.Cq, 04.50.-fh 

Recent observational data on type la supernovae explo- 
sions at high redshifts z = ^^^^ — 1 ~ 1 obtained indepen- 
dently by two groups [|l]j2| , as well as numerous previous 
arguments (see the recent reviews strongly support 

the existence of a new kind of matter in the Universe 
whose energy density not only is positive but also domi- 
nates the energy densities of all previously known forms 
of matter [here a{t) is the scale factor of the Friedmann- 
Robertson- Walker (FRW) isotropic cosmological model, 
and to is the present time]. This form of matter has a 
strongly negative pressure and remains unclustered at all 
scales where gravitational clustering of baryons and (non- 
baryonic) cold dark matter (CDM) is seen. Its gravity 
results in the present acceleration of the expansion of the 
Universe: a(to) > 0. In a first approximation, this kind 
of matter may be described by a constant A-term in the 
gravity equations as first introduced by Einstein. How- 
ever, a A-term could also be slowly varying with time. If 
so, this will be soon determined from observational data. 
In particular, if we use the simplest model of a variable A- 
term (also called quintessence in [|]) borrowed from the 
inflationary scenario of the early Universe, namely an ef- 
fective scalar field <& with some self-interaction potential 
U (<&) minimally coupled to gravity, then the functional 
form of ?7($) can be determined from observational cos- 
mological functions: either from the luminosity distance 
Dl{z) I^,^, or from the linear density perturbation in 
the dustlike component of matter in the Universe 5m{z) 
for a fixed comoving smoothing radius |6j. However, this 
model cannot account for any future observational data, 
in particular, for any functional form of Dl{z). This 
happens because a variable A-term in this model should 
satisfy the weak-energy condition e A +PA > 0. In terms of 
the observable quantity H{z) = a{t)/a{t) describing the 
evolution of the expanding Universe at recent epochs, the 



following inequality should be satisfied 

^^>3a„,oi?o'(i + ^f . (1) 

dz 

Here, Hq = H{z = 0) is the Hubble constant, rim,o is the 
present energy density of the dustlike (CDM-|-baryons) 
matter component in terms of the critical density £crit = 
3Hq/87tG {c = h = 1, and an index stands for 
the present value of the corresponding quantity). Note 
that the inequality (0) saturates when the A-term is ex- 
actly constant. It is not clear from the existing data 
whether (|l|) is satisfied at all. Actually the opposite 
holds: An attempt to reconstruct [/(<&) from the super- 
novae data and fitting of existing data to a model with 
a linear equation of state for the A-term p\ = we\, with 
w < —1 , shows that the possibility of violation of in- 
equality (|l|), though strongly restricted, is not completely 
excluded. Hence it is natural and important to consider 
a variable A-term in a more general class of scalar-tensor 
theories of gravity where the requirement (|l|) does not 
arise. Moreover, this generalization of general relativity 
(GR) is inspired by present more fundamental quantum 
theories, like M-theory. In these theories, the scalar field 
<f> is just the dilaton field, hence we shall call it so below. 

Thus, we are interested in a universe where gravity is 
described by a scalar-tensor theory, and we consider the 
Lagrangian density in the Jordan frame [l0| 

L = l{F{<i>) R - .g^^a^^a,*) - C/($) + i,„(.9^„) , 

(2) 

where Lm describes dustlike matter and F($) > 0. 
This corresponds to the Brans-Dicke parameter uj — 
F/{dF/d^Y > 0- may also introduce a function 



1 



Z($) in front of the kinetic term {d^^Y , but it can be set 
either to 1, or to —1 by a redefinition of the scalar field. 
Under the assumption of absence of ghosts in the theory, 
the second possibility requires the Brans-Dicke parameter 
to lie in the range —3/2 < a; < (see ||ll|] for more de- 
tails). Since this clearly contradicts solar system tests of 
GR either in the absence of U{^), or for U{^) satisfying 
the condition (||) below for scales of galaxies and clus- 
ters of galaxies, we will not discuss this possibility fur- 
ther. We do not introduce any direct coupling between $ 
and Lm (though this possibility could be envisaged, too). 
This guarantees that the weak equivalence principle is ex- 
actly satisfied (universality of free-fall of laboratory-size 
objects), and also that fundamental constants, like e.g. 
the fine-structure constant, do not change with time in 
this theory. This is in very good agreement with labora- 
tory, geophysical and cosmological data . 

Such a scalar-tensor theory was recently considered as 
a model for a variable A-term for some special choices 
of F{<^) and C/($) (see jl^). Our approach is just the 
opposite: We want to derive these functions from obser- 
vational data. Since we have to determine two functions 
F($) and U{^), we will need both observational func- 
tions Dl{z) and (5,„(z), in contrast to GR. Then the re- 
construction problem can be uniquely solved as will be 
shown below. Note that the angular diameter as a func- 
tion of z provides the same information as Dl{z) (see 
and the second reference in Q). 

It is most appropriate for us to work in the Jor- 
dan frame (JF), in which the various physical quanti- 
ties are those that are being measured in experiments, 
even though the Einstein frame (EF) often provides a 
better mathematical insight. In addition, the dilaton ap- 
pears to be directly coupled to dustlike matter in the EF 
frame, in contrast to the JF. For a flat FRW universe 
with ds^ = —dt^ -f a?dyp' , the background equations in 
the JF are then 



3i^ij2 = + _ + t/ _ 2,HF 
-2FH ^ pm + ^'^ + F - HF . 



(3) 
(4) 



Their consequence is the equation for the dilaton itself: 

$ + 3i/$+ ^ - 3(ij + 2i/2) ^ = . (5) 

Combining Eqs. (^)-(0) and changing the argument from 
time t to redshift z, we obtain the following basic equa- 
tion for F{z): 



F" 



{\nHy - 



1 + z 



F' 
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(1 + Z)2 
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2(lnH)' 

l + z 



F 



where the prime denotes the derivative with respect to z. 



The effective value of Newton's gravitational constant 
Gn in Eqs. (|-|) is given by the formula Gn = l/SnF. 
We shall use its present value Gn,q in the definition of 
the critical density Ecrit- On the other hand, Gn,o is not 
the quantity measured in laboratory Cavendish-type and 
solar-system experiments. For a massless dilaton, the 
effective gravitational constant between two test masses 
is given by 



Gcff — 



1 f2F + A{dF/d^f 



SttF \2F + 3{dF/d^)^ 



(7) 



In our case, the dilaton is massive, so the expression (|^) 
will be valid for physical scales R such that 



R ^ max 



d^U 



d^F 



d$2 



(8) 



Previously, the expression Goff was known from the post- 
Newtonian expansion; below we rederive it using the cos- 
mological perturbation theory. 

Let us now list the restrictions of the theory (^) which 
follow from solar-system and cosmological tests. The 
post-Newtonian parameters (3 and 7 for this theory are: 



7 = 1- 



{dFld<^f 



F + 2{dF/d^Y 



1 F {dF/d^) d-y 
^ ^ 4 2F iidF/d-^Y d$ 



(9) 
(10) 



Using the upper bounds on (7 
measurements |l^Jl^, we get 



1) from solar system 



Fo^ {dF/d^)l < 4 X IQ-^ 



(11) 



So, Gn^ and Gcft.o coincide with better than 2 x 10~^ 
accuracy. On the other hand, the difference between Gjq 
and Gcff may be larger at redshifts z 1 since neither 
the upper limit on /3, nor the present experimental bound 



IGefi/Gefil < 6 X 10^12 yj.-! ^ significantly restrict 



{d^F/d^'^)o- Note that we cannot use the nucleosynthe- 
sis bound on the change of Gcff since that time as the be- 
havior of Gcff during the intermediate period is unknown, 
unless we make additional assumptions (see below). 

The theory (^ describes a variable A-term with desired 
properties if the following three conditions are satisfied: 

1) The A-term is dynamically important at present, 
namely, f^A.o ~ 0.7 ~ 2Q.m,Q, or 



$2 

— + f/ - iHF 



0.7£c 



2p, 



m,0 



(12) 



2) The A-term has a sufficiently large negative pres- 
sure to provide acceleration of the present Universe. The 
condition ap > reads: 



2C/o > (Pm + 2$2 + 3^^ + iHF)c 



(13) 
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3) The dark matter described by the A-term remains 
unclustered at scales up to i? ~ 10h~^{l + z)~^ Mpc and 
probably even more (here h = TJq/IOO km s^^ Mpc^^). 
To achieve this, it is sufficient to assume that the in- 
equality (||) is satisfied for all scales in question. 

The first step of our program is purely kinematical: we 
determine H{z) from Dl{z) like in GR, 



1 



H{z) 



1 



(14) 



The functional dependence of Dl{z) on the cosmological 
parameters, like f2m,o, is of course model dependent. If 
rim.o is already known from other tests, we can find al- 
ready at that stage of the reconstruction a quantity such 
as the present effective equation of state of the dilaton 
from the formula (cf. ||]): 



Wo 



PA,0 



(2/3)(dlni7/dz)o - 1 

1 — flmO 



(15) 



e\^o contains the term —3HqFq, so that Om_o + ^A.o = 1- 
The dilaton equation of state can be determined for 
z > 0, too; one has only to define what should be called 
the pressure and the energy density of the dilaton in 
general. Actually, we will show below that ^m.o is it- 
self self-consistently determined from our approach, so 
no additional information is required to find w{z). 

In contrast to GR, Eq. (||) is no longer sufficient to de- 
termine U{z); one should know F{z), too. For this pur- 
pose we will use Sm{z). We consider perturbations in the 
longitudinal gauge ds^ = —(1 -I- 2(j))dt^ + a^(l - 2i/')dx^. 
Working in Fourier space (a spatial dependence exp(zk.x) 
with fc = |k| is assumed), the following equations are ob- 
tained: 



= i, = V - 5F/F , 

. e M + Hv) 

Om = ^ V + 3 

at 



(16) 
(17) 



where the gauge invariant quantity Sm = {5pm)/pm + 
3i?w, and v is the peculiar velocity potential of dustlike 
matter. We also get 

-ZF^ - {^^F - $2 + 2,hP^ 4) = 

F ■ f F^\ 
= PmS„, + 3j6F+i-- 6H^ - 3-^ j 5F 

+ ^S'<P + 3H<i> 6^ + SU, (18) 
and the equation for the dilaton fluctuations 5$: 



6^ + 3H5^ 



--3[H + 2H^) — 



d^F d^U 



(5$ = 



— (0 - 2i;) - 3(V^ + 'iHijj + H<j>) 



-{Sip + (^)$ - 20 



dU 



dF 
M 

(19) 



Let us now consider sufficiently small scales R = 2'Ka/k 
for which the inequality (^ is well satisfied. For exam- 
ple, if 5m{z) is determined from the abundance of rich 
clusters of galaxies, then the relevant comoving scale is 
R - 8/1^7(1 + z) Mpc. If the r.h.s. of Eq. (|) is ~ H^, 
then the corresponding small parameter is R^Hq ^ 10^^. 
Note that we have another parameter, uj^^, which is 
small at the present time, Eq. (pT|), but it need not be 
so small in the past. Also, this parameter may be larger 
than a^H'^/k'^. For this reason, we will first keep it. 



The solution of Eqs. (|T£-19) in the formal short- 
wavelength limit fc — > oo can be found following the an- 
alytical method used in [|j in the GR case, confirmed 
numerically in [Q. The idea is that the leading terms 
in Eqs. (|l6|~|l9|) are either those containing fc^, or those 
with 6m- Then, using (0) and the l.h.s. of Eq. (|l|), the 
standard form of the equation for dustlike matter density 
perturbation follows: 



2HS„ 



k 



2 -2 

a c 







(20) 



Now we consider the solution of Eq. ( p^ of interest to 
us, for which \5^\ ^ k'^a^'^\6^\. It corresponds to the 
growing adiabatic mode. So, keeping terms with fc^ in 
Eq. ( p^ and then using the r.h.s. of Eq. (|l6|), we obtain: 



5^ 



F dF/d^ 
F + 2{dF/d^y 



(21) 



In the GR case, 5$ oc k^'^cf) in the limit fc oo, so mat- 
ter producing the A-term is not gravitationally clustered 
at small scales (physically, due to free streaming). This 
is not so in scalar-tensor gravity: The dilaton remains 
partly clustered for arbitrarily small scales, this cluster- 
ing being small only because lu is large. 

Keeping only terms with fc^ or Sm in Eq. (|l^), we get 
the expression of 4> through Sm and SF. Finally, inserting 
it into Eq. (^o|) and using Eq. (|l|) , we arrive to the closed 
form of the equation for Sm '■ 

Sm + 2H5m - 47rGoff Pm Sm - , (22) 

with Gcff defined in (|^) above. In terms of z, (22) reads: 



3 



{l + z)H, 



2 ^cff(z) 



l + z 

f^m.O S. 



SL 



(23) 



Eq. (|2^) does not contain fc^ at all. Thus, its solutions, 
as well as the corresponding expressions for (5<I>, do not 
oscillate with the frequency k/a for fc — > oo. This justifies 
the assumption about (5$ made above. 

Extracting H{z) (from Di^{z)) and Sm{z) from ob- 
servations with sufficient accuracy, we can reconstruct 
Gcs{z)/Gn.o analytically. Since, as follows from 
Eq. (pi]), the quantities Gofr.o and Gn.q coincide with bet- 
ter than 0.02% accuracy, Eq. (23) taken at z = gives 



3 



also the value of flm,o with the same accuracy. Thus, 
in principle, no independent measurement of fim.o is re- 
quired. 

The resulting equation Gcff(z) = piz), where p(z) is 
a given function following from observational data, can 
be transformed into a nonlinear second order differential 
equation for F{z) if we exclude d$ (which appears in 
dF/d^) using the background equation (|^), which reads 
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$'2 = 



(IniJ)' 



z 



F' 



2{\tlH)' 

l + z 



F-3(l 



m.O 



(24) 



Therefore, F{z) can be determined by solving that equa- 
tion provided Fq (— I/SttGna) and F^ are known. 

However, this procedure can be simplified a lot un- 
der reasonable assumptions, and taking into account 
the small present values of ui^^ = F^^{dF/d^)^ and 
Geff/Geff. Indeed, the value of uj^^ for < z ^ 1 can 
be estimated from the first terms of its Taylor expansion 
lOq^ + z {duj~^ / dz)Q. Neglecting contributions propor- 



tional to ujq , we then get u;-^ ~ 2 z A {d^F/d^^ 



with 



A = -(dlnF/d$)o$o/i?o, whereas Goff/Geff ~ AiJo [1- 
{d^ F / d^'^)Q\. If {d'^F/d<i>^)o diflPers significantly from 1, 
we can thus conclude that < |2Gcfr/-ffoGoff| < 0.25. 
On the other hand, if {d^F/d^^)o happens to be close to 
1, one can still assume that there is no special cancella- 
tion of large terms in the r.h.s. of Eq. (^), and therefore 
that ^ GFqHq. The above estimate for u;^^ then gives [g] 
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